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With  the  increasing  complexity  and  astronomical  costs 
associated  with  tho  dosign  and  production  of  spaco  systems, 
availability  and  reliability  are  among  the  most  important 
factors  to  consider  when  building  space  system*.  Also,  as 
repair  of  space  systems  becomes  more  feasible  with  the  advent 
of  the  shuttle  and  other  upcoming  space  projects, 
maintainability  is  becoming  another  important  factor  in  space 
system  design.  But  none  of  these  factors  can  be  applied 
without  a  thorough  understanding  of  the  theory  behind 
availability,  reliability,  and  maintainability. 

The  purpose  of  this  study  is  to  verify  the  hypothesis 
that  Monte  Carlo  simulation  is  the  best  method  of  finding 
lower  confidence  bounds  for  the  availability  and  reliability 
of  maintained  systems  in  different  configurations.  Also, 
this  study  verifies  that  once  a  lower  confidence  bound  is 
found,  it  is  the  same  over  the  lifetime  of  a  maintained 
system . 

I  would  like  to  thank  cay  thesis  advisor,  Professor 
Albert  K.  Moore,  for  his  most  valuable  and  very  expert  advice 
and  guidance  during  this  study.  I  would  also  like  to  thank 
Lt  Col  Joseph  W.  Coleman,  my  reader,  for  hi*  help  and 
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guidance  during  this  study-  Finally,  I  wish  to  nhank  ray 
husband,  Nick,  for  his  encouragement  and  moral  support  during 
this  study  and  throughout  this  master's  program- 


Karan  S.  harland 
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This  research  determined  the  feasibility  and  sfficisncy 
of  a  Monts  Carlo  method  of  simulating  the  lower  confidence 
limits  for  the  availabilities  and  reliabilities  of  maintained 
systems.  The  steady-state  availabilities  of  single-unit 
systems  and  the  time-constrained  availabilities  and 
reliabilities  of  two-unit  parallel  systems  were  simulated. 

Q—  — ' 

First,  a  baseline  of  "true*  exponentially-distributed 
Mean  Time  Between  Failures  (MTBFs)  and  Mean  Time  To  Repairs 
( MTTRs )  were  simulated  using  the  chi-square  distribution. 

Then  other  MTBFs  and  MTTRs  were  simulated  to  represent 
sampling  of  other  systems.  The  availabilities  and 
reliabilities  were  found  using  these  simulated  MTBFs  and 
MTTRs.  Next,  simulated  availabilities  and  reliabilities  were 
ordered,  and  lower  confidence  limits  v?ere  found.  These  lower 
confidence  limit  point  estimates  were  compared  against  the 
systems '  exact  availabilities  and  reliabilities.  Lastly,  the 
success  of  this  Monte  Carlo  method  is  determined  by  how  well 
the  simulated  lower  confidence  limit  availability  and 
reliability  point  estimates  cover  the  exact  availabilities 
and  reliabilities. 


xx 


VERIFICATION  07  A  MONTE  CARLO  SIMULATION 
METHOD  TO  FIND  LOWER  CONFIDENCE  LIMITS 
FOR  THE  AVAILABILITY  AND  RELIABILITY  OF 
MAINTAINED  SYSTEMS 

I.  Introduction 

The  reliability  and  maintainability  disciplines  became 
firmly  established  in  the  1950's  due  to  the  United  States' 
increasing  emphasis  on  its  military  and  space  programs.  The 
need  for  engineers  to  develop  more  efficient  and 
longer-lasting  equipment  for  these  programs  was  due  to 
several  f actors t 

a)  high  equipment  failure  rates 

b)  inflation  which  increased  the  cost  of  buying  and 
maintaining  equipment 

c)  increased  equipment  complexity 

d)  the  desire  to  develop  a  methodical  approach  to 
minimize  causes  of  failure  (Fullerton,  1969:1). 

Then  and  now,  engineers  are  interested  in  three  general 
areas  when  designing  systems  -  reliability,  maintainability, 
and  availability.  Reliability  is  "the  probability  that,  when 
operating  under  stated  environmental  conditions,  the  system 
will  perform  its  function  adequately  for  a  specified  interval 
of  time"  (Kapur ,Lamberson,  1977:1).  Maintainability  is  "the 
probability  that  a  failed  system  can  be  made  operable  in  a 
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specified  interval  of  downtime"  ( Kapur  ,  Laoberson ,  1977:225). 
Thus,  a  maintained  system  is  one  that  can  be  repaired.  The 
measure  that  includes  both  reliability  and  maintainability  is 
availability.  Availability  is  "the  probability  that  a  system 
is  operating  satisfactorily  at  any  point  in  time  and 
considers  only  operating  time  and  downtime,  thus  excluding 
idle  time"  < Kapur , Lamberson ,  1977 i 225).  An  engineer  can  be 
sure  his  design  is  sound  if  he  can  measure  the  reliability 
and  availability  of  the  system,  and  those  measures  fall  into 
his  specified  confidence  interval.  A  confidence  i  '•erval  is 

a  range  of  values  which  is  believed,  with  >t 
preassigned  degree  of  confidence,  to  include  the 
particular  value  of  some  parameter  or  charac¬ 
teristic  being  estimated.  The  degree  of  con¬ 
fidence  is  related  to  the  probability  of  ob¬ 
taining  by  random  samples  ranges  which  are 
correct  [ James , Beckenbach ,  1968:65]. 

Many  statisticians  have  done  previous  work  in  developing 
and  evaluating  techniques  to  obtain  the  confidence  limit*  of 
system  reliabilities  and  availabilities  modeled  by  different 
underlying  distributions.  One  of  the  first  works  presented 
came  from  Mary  Thompson.  Thompson  developed  analytical 
techniques  to  determine  confidence  limits  for  the 
availability  of  exponentially  distributed  repairable  systems 
in  1966  (Thompson,  1966).  In  the  following  year,  Louis  Levy 
nd  Albert  Moore  introduced  a  Monte  Carlo  simulation 
technique  to  estimate  reliability  confidence  limits  from 
component  test  data  with  normal,  gamma,  or  Weibull 


probability  distributions  (Levy, Moore,  1367).  H.  L.  Gray  and 
W.  R.  Schucaay  extended  Thompson's  techniques  to  find 
availability  confidence  limits  for  systems  with  lognormally 
distributed  repair  time*  in  1969  (Gray ,Schucany ,  1969). 
However,  work  has  >t  to  be  done  to  obtain  confidence  limits 
fo  :he  availability  and  reliability  of  maintained  systems 
usi^g  Monte  Carlo  simulation  with  simulated  system  data. 

Statement  of  the  Problem 

Markov  models  and  exact  analytical  techniques  have  been 
used  to  produce  confidence  limits  for  the  reliability  and 
availability  of  maintained  systems.  However,  an  system  size 
and  complexity  grows,  these  models  and  techniques  become  very 
hard  to  apply.  A  quick  and  efficient  Monte  Carlo  simulation 
method  is  needed  to  find  confidence  limits  for  the 
reliability  and  availability  of  maintained  systems,  and  an 
assessment-  of  the  accuracy  of  the  method  is  needed.  According 
to  current  literature,  a  Monte  Carlo  simulation  can  identify 
independent  or  dependent  system  failures  occurring  in  either 
simple  or  complex  systems  (Almassy,  1979:366).  .iont®  Carlo 
simulation  can  also  model  systems  with  general  failure  and 
repair  distributions.  (In  this  case,  the  failure  and  repair 
distributions  are  assumed  to  be  exponential  because  the 
exponential  distribution  accurately  models  the  failure  of 
most  complex  electronic  ay-uterus  and  provides  a  fair  uodel  f  or 
the  repair  of  electronic  systems  [Thompson ,  1966:363- ) 


Therefor*,  Mont*  Carlo  simulation  will  provide  an  ®a»y  and 
accurate  method  of  obtaining  reliability  and  availability 
confidence  limits  for  maintained  systems. 


flfeiSS&LZim  fit  till:  R«g«ftifgh 

The  overall  objective  of  this  research  is  to  develop  ax 
efficient  Monte  Carlo  simulation  to  find  confidence  limits 
for  the  reliability  and  availability  of  maintained  systems. 
The  subobjectives  for  this  research  are  as  follows: 

a)  to  develop  a  Mont*  Carlo  simulation  method  to  find 
availability  confidence  limits  for  small  sample, 
single  system  configurations 

b)  to  find  the  availability  for  a  single  system  con¬ 
figuration  system  configuration  using  an  exact 
analytical  technique 

c)  to  evaluate  the  accuracy  of  the  Monte  Carlo 
simulation  method  by  comparing  the  Mont*  Carlo 
confidence  limit*  to  the  exact  availability 
for  <a  single  system  conf  igvirat  i  on 

d)  to  expand  the  Mont#  Carlo  simulation  to  get  both 
reliability  and  -availability  confidence  limit*  for 
small  sample,  parallel  system  conf igurat ion*  wioh 
time  and,  if  time  permit*,  more  complex  systems 

e  >  to  present  t,  «  actual  percentage  covsrag*  of  tha 
confidence  Internals  of  the  complex  systems  as 
another  wav  of  determining  the  accuracy  of  +l'i* 

•% 


method  of  obtaining  confidence  limits. 

Mt&lSdSlaSlX 

The  overall  approach  for  this  research  is  to  use  a  Mont® 
Carlo  simulation  to  estimate  the  availabilities  and  relia¬ 
bilities  of  maintained  systems  and  determine  ths  lower  confi¬ 
dence  limits  fc  oth.  This  Monts  Carlo  simulation  will 
generate  samp Is  Mean-Time-To-Repair  (MTTR)  and  Mean-Time-Be- 
tween-Fail  \xm  (MTBF)  values  from  the  chi-squar®  distribution. 
The  MTTRs  .  nd  MTBF*  are  exponentially  distributed  and  arc 
from  small  >  amplss  <t*n  to  fifty  systsm  configurations). 

Ths  saaplss  sr*  Type  II  censored  samples  which  are  a  set  of 
n  items  that  are  tested  until  a  specific  number  of  failures 
have  occurred  (Moore,  1963s  460).  Input  sample  MTTR  (  <6  )  and 
MTBF  (  9  )  values  will  form  the  "real  world"  baseline  data 
against  which  the  rest  of  the  values  will  be  compared. 

The  simulation  will  initially  consist  of  two  parts. 
First,  only  the  availability  of  a  single  system  configuration 
will  be  estimated  using  the  generated  MTTRs  and  MTBFs.  The 
estimated  availability  is  found  from  the  following  equation 
( Kapur , Lamberson ,  1977;2ZB) 

As»  "  MTBF/ (MTBF  +  MTTR)  (1) 

The  reliability  of  a  single  system  configuration  will 
not  be  estimated  because  if  a  system  is  not  available,  it  has 
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failed,  and  at  that  point,  it  ha*  reached  its  reliability 
limit.  However,  for  the  second  part  of  the  simulation,  both 
the  availability  and  reliability  of  two-unit  parallel  systems 
with  repair  will  be  estimated,  again  by  using  the  generated 
MTTR*  and  MTBFs.  The  availability  and  reliability  equation* 
for  parallel  systems  am  as  follows  tShooman,  1966 s 341-346 > : 

A<t>  -  <l-<AA/r3r4) )  -  XA^r3r4  <2) 

,_r  t,_  r  t ,  „  . 

3  /r3  “  #  4  /r4> 

where 

A  -  1/0 

Li  »  1/0 

)i~  2  A  for  an  ordinary  system 
A7*  A  for  a  standby  system 
fj  »  n  for  one  repairman 
t*  m  for  more  than  one  repairman 

<k1  >  1) 

n 

H  »  u  for  one  repairman 
U  *  7. li  for  two  r«p«irn«a 


•f  . 

M  f or  norm  than  two  repairman 

(k2  >  2) 

r3'  r4  “  £"<  X  +A/+^/+^v)  +/- 


(  ( <  X  +  A'-iV*  U)**2) 

-  <4t  AA'+A'jA  AfV>  )1/Z]/2 


where 


t(t)  *  (  A  +  ^  +  r^/r^-r^)**^ 

--  <  A  +  //+  r2/rl“r2 ***2^ 

+  <*''vvV 

-  (  A'/r1-r2)er2t 


t,r2  -  C-<>+A/+^/>  ♦/“  ( (  A  +  A'+  n')  **2 

-  4AA'  )1/2J/2 

A  ■  2  A  for  an  ordinary  system 

A  »  A  for  a  standby  system 

/j' m  I*  for  one  repairman 
/ 

M  *  k;i  for  mere  than  one  repairman  (k  >  1 


Reliabilities  can  be  found  for  the  two-unit  system 
configuration  because  one  unit  In  a  parallel  configuration 
may  fail,  but  ite  parallel  unit  may  still  be  available,  so 
bhe  configuration  still  works.  Thus,  the  parallel  system 
configuration  has  a  reliability  that  differs  from  the 
availability.  <The  repairable,  parallel,  single  repairman 
configuration  will  be  considered  in  this  simulation;  others 
will  be  considered  as  time  permits.) 
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After  both  parts  of  the  simulation  have  been  run,  the 
lower  confidence  limits  are  found  and,  by  simulating  a  large 
number  of  trials,  a  higher  level  of  accuracy  in  estimation 
can  bm  obtained.  This  Monte  Carlo  simulation  method  can  be 
assessed  and  verified  by  how  times  the  confidence  interval 
contains  the  true  system  availabilities  anc  reliabilities. 

Sumtaai  al  EsaMftfci&lflft 

The  next  chapter  summarizes  the  literature  available  on 
this  subject  from  the  1960s  to  the  present.  Chapter  3 
presents  the  Monte  Carlo  simulation  models  designs  and  the 
equations  and  variables  used  in  the  model.  Chapter  4 
discusses  the  simulation  rssults  in  answer  to  the  thesis 
subobjectives  and  conclusions.  Lastly,  Chapter  5  presents  a 
summary  of  the  research  effort  and  recommendations  for  future 
studies . 
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II-  ki&fc£A&L££ 

As  previously  stated  in  Chapter  I ,  statisti clans  and 
theorists  have  been  developing  various  techniques  to  find 
availabilities  and  reliabilities  of  many  different  system 
configurations  with  many  different  underlying  failure  or 
repair  distributions.  This  chapter  will  review  some  of  that 
work,  divided  into  these  different  areas  as  such  as  possible! 
availability,  reliability,  Monte  Carlo  simulation  of 
availability  or  raliability,  and  classic  books. 

H.asJsji  £&  AayjAfcillSJL 

There  have  been  several  works  published  focusing 
stzictly  on  availability  since  the  mid-1960s.  In  the  year 
following  Thompson's  work  on  analytically  determined 
availability  confidence  limits,  John  Buzacotu  showed  how  to 
reduce  reliability  block  diagrams  of  repairable 
series-parallel  systems  to  find  the  MTBF  and  availability  of 
the  system  components  (Buzacott,  1967).  Also  in  1967,  H.  L. 
Gray  and  T.  0.  Lewis  presented  an  analytical  method  to  find 
exact  system  availability  confidence  limits  with 
**l onentially  distributed  MTBFs  and  lognormrlly  distributed 
MTTRs  (Gray, Lewis,  1967).  In  1969,  Kenneth  Grace  gave  a 
comparison  between  a  Markov  model  and  some  approximation 
techniques  used  to  find  the  steady--ztate  availability  of 
repairable  systems  with  limited  component  spares  (Graco, 


1969).  In  1970,  Buxacott  did  some  mor®  r®»«arch  on 
repairable  systems,  this  time  using  nstwcrk  diagrams  to  find 
system  availability  and  failur®  frequencies  (Buxacott, 

1970b).  In  1973,  David  Chow  presented  a  math®mat iaal  modal  to 
find  tha  availabilities  of  a  redundant,  repairable  system 
with  a  standby.  The  system  had  to  have  q  components 
available  to  work,  and  the  underlying  failure  rates  were 
constant  with  generally  distributed  MTTRs  (Chow,  1973). 
Several  years  later  and  building  on  the  works  of  Thompson, 
Gray  and  Lewis,  and  Gray  and  Schuaany,  Mohamad  Hasaballa, 
Albert  Moore,  and  Joseph  Cain  presented  an  exact  analytical 
method  of  finding  lower  confidence  limits  on  steady-state 
asymptotic  availabilities  of  systems  with  exponentially 
distributed  MTBFs,  and  exponentially  and  lognormally 
distributed  MTTRs  (Hasaballa  and  others,  1983).  In  1983, 
another  work  concerning  repairable  systems  was  published  by 
Richard  Kenyon  and  Richard  Newell.  They  gave  an  exact 
analytical  solution  and  a  Fortran  program  designed  to  t'-£ 
the  steady-state  availability  of  a  system  needing  k-out-of-n 
components  to  function,  but  allowing  only  one  repair  fox  the 
systom  lifetime  ( Kenyon , Newell ,  1583).  During  the  same  year, 
JEtpaii  Funnemark  and  Bent  Matvig  broke  new  ground  with  th®i  r 
presentation  on  how  to  find  upp®r  and  lower  confidence  limits 
for  the  availability  of  multi-state  component  system* 

(  Funnemark  ,  Ma  tv  i  q  ,  1983).  Lastly,  in  l'3J4t  Ignas  10  M«nd*s 
introduced  a  ®»thod  of  finding  tn#  steady-state  availability 
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for  n  unlike  parallel  or  n  unlike  nerias  component*  with 
exponentially  distributed  failure  and  repair  rates  (Mender, 
1384).  All  the  above  works  have  bearing  on  this  thesis 
because  they've  considered  repairable  or  complex  systems, 
exponentially  distributed  failure  or  repair  rates,  confidence 
limits,  or  steady-state  availability.  The  next  section 
presents  reliability  works  which  have  considered  the  same 
f  actors . 

Maeltl  £&  Rsli ability 

The  list  of  theorists  working  on  reliability  prediction 
techniques  is  quite  extensive,  and  this  section  gives  only 
the  highlights.  Starting  in  the  early  1360s,  Thomas  Burnett 
I©  and  Beverly  Wales  presented  an  analytical  method  and  a  simple 

Monte  Carlo  method  to  determine  system  reliability  confidence 
limits  from  component  test  results  ( Burnett , Wales ,  1961).  In 
1963,  Oscar  Bernhoff  researched  a  method  bo  analytically 
determine  system  reliability  confidence  limits  by  combining 
the  components'  reliability  estimates  to  get  the  overall 
failure  probability  density  function.  Ho  also  determined 
that  a  Monte  Carlo  simulation  of  the  overall  system 
distribution  was  much  easier  to  use  when  component  estimates 
were  from  two  or  more  dissimilar  distributions  (Bernhoff, 
1963).  Also  in  1963,  Malcolm  McGregor  worked  out  a  method  of 
determining  the  reliabilities  of  repairable  systems  with  n 
identical  parallel  components  having  exponentially 
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distributed  MTBFs  and  MTTRs  (McGregor,  1963)  „  In  1963,  Albert 
Madansky  revealed  how  to  combine  component  reliability 
estimates  gathered  by  separate  tests  to  determine  the  overall 
complex  systam  reliability  (Madansky,  1963).  Then,  in  1966, 

L.  Mtun  used  McGregor's  findings  in  his  work  about  using 
transition  diagrams  describing  different  states  of  repairable 
or  nonrepairable  systems  to  determine  reliability  (Htun, 
1366).  In  1970,  Buiacsott  presented  yet  .mother  study,  this 
time  on  special  Markov  techniques  for  determining  the 
reliability  (and  availability)  for  a  large  number  of  states 
in  repairable  systems  (Buxacott,  1970a).  (Earlier,  in  1966, 
Kenneth  Blakney  and  Frederick  Dietrich  presented  a  thesis 
covering  Markov  reliability  processes  known  to  that  time;  R. 
Fullerton  updated  thit  list  in  1969  [Blakney , Dietrich ,  1966; 
Fullerton,  1D69].)  Then,  in  1972,  building  on  Madansky's 
work,  Robert  Easterling  presented  a  technique  for  determining 
the  maximum  likelihood  estimate  of  system  reliability  from 
component  test  results  and  substituting  that  estimate  into  an 
incomplete  bet*  funwtio"  to  determine  reliability  confidence 
limits  (Easterling,  19 n7 t  ,  Next,  Nancy  Mann  and  Frank  Grubbs 
presented  several  methods  of  approximating  reliability 
confidence  limits  for  series  or  parallel  systems  using 
component  test  results  gathered  from  Type  I  or  Type  II 
censored  sampling;  they  based  their  work  on  studies  done 
previously  by  Madansky  and  Easterling  ( Mann , Grubbs ,  1974). 

Later,  in  1977  and  1978,  Edward  Bil’kam  and  Albert  Moore 


published  a  extended  study  on  estimating  reliability 
confidence  limits  from  multiple  independent  grouped  censored 
samples  with  failure  times  known  or  unknown  (Bilikam, Moore, 
1977  and  1978).  Looking  again  at  the  technique  of  using 
Markov  processes  to  determine  system  reliability,  Joseph 
Foster  and  Alberto  Garcia-Diaz  formulated  generalized  Markov 
models  for  three  classes  of  reliability:  systems  with 
catastrophic  failure,  systems  that  must  be  down  to  be 
repaired,  and  systems  that  can  be  repaired  while  functioning 
(Foster, Garcia-Diaz,  1982).  Again,  in  the  area  of 
multi-state  systems,  David  Butler  presented  mathematical 
computations  for  finding  complex  multi-state  system 
reliability  bounds  in  1982  (Butler,  1982).  In  the  same  year, 
I.  Gertsbakh  also  presented  mathematical  formulas  to 
determine  upper  reliability  confidence  limits  for  parallel, 
series-parallel,  and  k-out-of-n  systems,  all  with  exponential 
component  distributions,  and  Type  I  or  Type  II  censored 
sampling  (Gertsbakh,  1982).  More  mathematical  formulas  were 
presented  by  Tetsuo  Miyamura  for  combining  exponentially 
distributed  component  and  system  estimates  to  determine 
component  reliability  (Miyamura,  1982).  Lastly,  one  of  the 
more  recent  works  was  presented  by  Tze  Li  in  1984  when  he 
demonstrated  the  use  of  an  empirical  Bayes  estimator  for 
determining  system  reliability  using  a  large  sample  size  of 
systems  with  exponentially  distributed  failure  rates  (Li, 

198  4  )  . 
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Monte  Carlo  Computer  Simulation  Works 

The  projects  presented  in  this  section  represent  work, 
done  in  the  areas  of  availability  or  reliability  specifically 
using  Monte  Carlo  simulation  techniques.  Beginning  in  1960, 
Donald  Orkland  published  a  pioneer  simulation  study  on 
finding  lower  confidence  limits  for  system  reliability  using 
only  sample  component  failure  data  (Orkland,  I960).  Two 
years  later,  W.  Connor  and  W.  Wells  expanded  Orkland #s  work 
with  a  study  of  their  own  which  found  system  reliability 
confidence  intervals  for  serially  arranged  components  with 
binomially  distributed  failure  rates  ( Connor , Wells ,  1962). 
Next,  in  1964,  Louis  Levy  devised  a  Monte  Carlo  method  to 
find  system  reliability  confidence  limits  using  component 
test  data  with  exponential,  WeibulX,  gamma,  normal,  or 
lognormal  failure  rates.  He  based  his  method  on  previous 
works  by  Burnett  and  Wales,  and  Bernhoff  (Levy,  1964).  The 
following  year,  Albert  Moore  presented  an  extension  of  the 
Monte  Carlo  technique  used  by  Levy  to  find  confidence  limits 
when  the  distribution  or  joint  distribution  of  the  estimators 
were  unknown  and  the  data  came  from  Type  I  or  Type  II 
censored  sampling  (Moore,  1963),  Then,  in  1967,  Levy  and 
Moore  presented  a  joint  paper  on  a  more  efficient  Monte  Carlo 
simulation  technique  to  estimate  reliability  confidence 
limits  from  component  test  data  having  normal,  gamma,  or 
Weibull  failure  distributions  (Levy, Moore.  1967).  In  1968, 
Leonard  Doyon  and  Martha  Berssenbrugge  presented  a 
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computational  method  which  represented  system  states  as 
differential  equations  and  iteratively  calculated  reliability 
and  availability  estimates;  they  based  their  work  partly  on 
Htun's  work  (Doyon , Bersssnbrugge ,  I960).  During  that  same 
year,  Donald  Gilmore  presented  a  study  that  used  Monte  Carlo 
methods  to  simulate  failures  of  a  component  reliability  block 
diagram  to  estimate  the  overall  complex  system  reliability 
(Gilmore,  1968).  In  1972,  Robert  Lannon  devised  a  Monte 
Carlo  technique  to  determine  the  reliability  of  a  complex 
system  with  dissimilar  components  having  Weibull  failure 
distributions.  Previous  studies  by  Bernhoff,  and  Levy  and 
Moore  provided  a  basis  for  his  work  (Lannon,  1972),  Next, 
over  the  period  of  two  years,  Satish  Kamat,  Michael  Riley, 
and  William  Franzmeier  jointly  published  studies  on  using 
Monte  Carlo  simulation  to  find  complex  system  reliability  by 
inputting  reliability  flaw  graphs  and  minimal  tie-sets 
( Kamat , Riley ,  1975;  Kamat , Franzmeier ,  1976).  Next,  Hiromitsu 
Kumamoto,  Kazuo  Tanaka,  and  Koichi  Inoue  published  a  study 
based  on  previous  work  done  by  Levy  and  Moora,  and  Kamat, 
Riley,,  and  Franzmeier.  They  developed  a  Monte  Carlo 
simulation  to  find  system  reliability  using  reliability  block 
diagrams  or  fault  tree  analysis  with  a  controlled  variate  or 
importance  sampling  technique  (Kumamoto  and  others,  1977).  In 
1980,  Malcolm  Easton  and  C.  Wong  used  a  sequential 
destruction  method  to  est’mate  reliability  for  systems  of  100 
or  more  components  and  wi\n  dependent  or  independent 


component  failures;  their  work  is  based  on  the  work  of 
Kumamoto,  Tanaka,  and  Inoue  ( Easton , Wong ,  1980).  In  turn  and 
in  the  same  year,  Kumamoto,  Tanaka,  Inoue,  and  Ernest  Henley 
presented  a  study  basQd  on  Easton  and  Wong's  work  which  found 
the  unavailability  of  large,  complex  systems  by  using  dagger 
sampling  which  decreased  the  number  of  trials  and  computation 
time  of  the  Monte  Carlo  simulation  'Kumamoto  and  others, 
1980a).  They  also  published  another  study  in  late  1980  which 
gave  a  technique  to  estimate  the  unreliability  of  large, 
repairable  systems  using  a  state  transition  or  indirect  Monte 
Carlo  method  (Kumamoto  and  others,  1980b).  Two  years  later, 
another  foursome  presented  a  study  on  three  different  Monte 
Carlo  techniques  used  to  find  system  reliability  confidence 
limits;  Kathleen  Depuy ,  Jon  Hobbs,  Albert  Moore,  and  J. 
Johnston  developed  and  analyzed  the  accuracy  of  univariate 
asymptotic,  bivariate  asymptotic,  and  'modified  double  Monte 
Carlo'  techniques.  They  based  their  work  o*?  studies  done  by 
Orkiand,  and  Levy  and  Moore  (Depuy  and  others,  1982).  Tn 
1983,  Roy  Rice  and  Albert  Moore  published  a  paper  describing 
a  Monte  Carlo  simulation  using  pass-fail  data,  SDd  they 
explained  how  to  treat  no  failure  cases.  Their  work  was 
baaed,  in  part,  on  Easterling's,  Orkland's,  and  Levy  and 
Moore's  previous  works  (Rice,Moor«.  1983).  Lastly,  one  of 
the  latest  works  was  presented  in  April,  1385  by  Kadaba 
Gopalakr l shnan  who  developed  a  Monte  Carlo  model  to  find 
general  system  reliability,  maintainability,  and  availability 


IB 


with  interruption  factor*  and  delays  in  repair  taken  into 
account.  His  model  used  findings  previously  given  by 
Gilmore,  Kamat, Riley  and  Franrmeier,  and  Buxacott 
(Gopaiakrishnan,  1985). 
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CllMAsaJL  Bsok* 

Although  not  previously  mentioned,  the  ba«,i£  for  .aany  of 
the  articlw*  on  availability,  reliability,  and  Wont#  Carlo 
simulation  summarised  above  cam*  from  two  hooks.  The  first 


bo*, 4r. , 


itv .  was  written  by 


Richard  Barlow  and  Frank  Proschan  in  1965  and  presented  basic 
mathematical  theories  for  reliability  ( Barlow , Proschan , 


1965).  The  second  book. 


: _ An. 


Engineering  Approach-  was  written  by  M.  L.  Shooman  in  1968 
which  gave  an  engineering  perspective  to  reliability  (i.e., 
now  to  apply  reliability  in  engineering)  (Shooman,  1968). 
Both  books  are  classics  because  they  helped  lay  the  basic 
groundwork  from  which  other  studies  were  generated. 


Conclusion 

Many  studies  have  been  done  in  the  areas  of  system  and 
component  availability  and  reliability,  but  many  more  are 
<r#  to  be  don.  as  today's  equipment  and  systems  grow  even 
max  e  complex  . 
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This.  chapter  d«?arib*»*  th*  d®«ign  and  coaput&tion®  u»«d 
in  th*  Won't*  r*rlo  sisulation  aod«l .  Th«  simulation  mod*l 
wan  u*«d  to  find  th®  *va&iiabiliti®«  of  *lngl«-unit  systems, 
and  th®  *vailabiliti*s  and  rsliabilit  ;.®s  of  two-unit  ordinary 
parallel  syst •*«. 


'  i 


c 


ilmli  SxMsa 


Computations .  Th®  simplest  *jy*t«*  to  analys®  is  a 
single-unit  syst®a,  This  Mont®  Carlo  model  simulates  th® 
availability  <A(t}  when  t  is  tim®>  of  singl®~uni t  systems 
havi.ng  an  exponential,  or  constant,  failur®  distribution®. 
Now,  if  th®  single-unit  sys ifc«sa  ha*  no  capability  to  be 
repaired,  the  system  availability  will  aqua!  th®  ®y*t«i«i 
r®l  iability .  But,  having  a  rapaii;  capability  incraases  th® 
system  availability  -  *  desirable  syst«t&  character i»ti o.  Th® 
mod®l  computation  calculates  th®  availability  for  a 
repairable  system  which  will  bn  l»!'g«r  them.  th®  ®y*t»m 
reliability.  A  Martov  graph  mod®l*  th®  states  of  * i ngl«-~uni t 
system  availabilities  in  Figure  3.1  (Sh.  oraan,  136S<  33B) . 

Th®  symbols  in  Figur®  3.1  ar* 

* i  »  tn*  (or  systoa}  in  a  a tat® 

A  33  1  /M®an  Tim®  Between  Failure  (HTL'F) 
fi  *  1/Mean  Tim®  To  R®j>a.ir  ■  MTTR) 


r 


Ifi 


i»"  p 


0 


■  the  operating  state  of  the  system 
•1  ■  the  failura  state  of  tha  systaa. 

A  t  *  change  in  tiaa 

(If  tha  aystaa  was  not  repairable,  s^  would  ba  an 
“absorbing"  state ,  and  no  return  to  Sq  would  be  possible  .) 


Figure  3.1  Markov  Graph  for  Single-Unit  Systaa 
Availability  (Shooman,  1968 i 338) 
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The  Markov  gr»ph  sbttte  probability  •^w  at  ion*  used  in 


this  si»ulaticm  Model  Mare  a*  follows; 

p®0Ct)  *  it/ik  +  n)  +  <  k /{  k  +  V  Ms*p<-(  X  +  M  Ht)  > )  >  <4> 

F*1<t)  «  X/ <  X  +  ju  >  -  (  X  /<  X  +  w  )  l«p( -<  X  +  fi  >  (t) ) ) )  (5) 

So,  by  definition,  the  availability  (or  probability  of  the 
system  being  operational)  is  PsQ(t). 

A(t)  *  PsQ(t)  (6) 

Mow,  ass  t  (or  time)  gets  larger,  the  availability 
function  approaches  approaches  a.  certain  value  -  a 
“ steady-state*’  value.  For  the  single -unit  system,  the 
steady -state  availability  (A  (t))  is 

3S 

As*u)  -  lia  A(t)  -  /</(  X  +  )  (7) 

Equation  7  is  the  equation,  used  in  the  Monte  Carlo  model 
to  compute  the  steady -state  availabilities  of  sample 
single-unit  systems  (Shoo  <ian,  1968:338). 


Model .  Th®  simulation  modal  used  th®  following  baaalin® 

Mean  Tia®  B®tw«®n  Failur®  (1/  X)  -  100  hour* 

Mean  Tim®  To  Repair  (1/ (J  )  »  20  hour* 

Fir*t,  n®w  ®®timat®s  of  MTBF»  (9)  !»nd  MTTR*  (  (f>  )  wer® 
g®n®rat®d,  thus  stimulating  MTBFe  and  MTTfta  from  "tru®" 
aingl®— unit  *yst®m*.  Th®*®  new  MTBFs  and  MTTR*  were  th®n 
u*®d  to  simulat®  new  sampl®  MTBFs  and  MTTRs.  Both  *ampl®a 
were  drawn  using  the  chi-square  distribution  as  follow*: 


8  *  <xZ2r9>/2r  (8) 

where 

r 

«  number  of  failures  for  th®  Typ®  II 

c«n*ored  sample 

X2 

X  2r 

*  random  number  drawn  from  the 

chi-square  distribution  with  2r  degrees 

of  freedom 

and  MTTRs  were  generated  using  the  same  equation  with 
substituted  for  0. 

Then ,  each  sample  6  and  <f>  were  converted  to  A  and  // 
respectively  (  A  ~  1/0  and  //  *  1/0).  Finally,  the 
simulated  steady-state  availability  point  estimates  were 


found  u»ing 


Aaa<t>  «  n/{  K +  U  )  ( 9 ) 

In  summary,  simulated  steady-state  availabilities  ware 
found  for  Type  II  censored  sample  sizes  of  10,  20,  30,  and 
SO.  Point  estimate*  war*  found  for  the  low*r  confidanaa 
limit*  of  0.98,  0.9S,  0.90,  0.8S,  and  0.80.  Tha  modal  ran 
for  S00  rapatitions  with  300  trials  par  rapatition. 

Two-Unit  Parallal  System 

This  Monta  Carlo  modal  simulates  tha  availabilitias 
( A  ( t ) )  and  reliabilities  (R(t) )  of  samples  of  two-unit 
ordinary  identical  parallel  systems.  The  system*  have  both 
exponential,  or  constant,  failure  and  repair  rates.  The 
computations  and  modal  for  both  availability  and  reliability 
are  discussed  below. 

Availability  Computation .  Due  to  tha  nature  of  it's 
configuration,  the  parallel  system  can  have  a  repair 
capability  that  increases  it's  overall  availability.  If  one 
unit  in  tha  system  fails,  tha  system  can  keep  functioning 
with  the  operational  unit  while  the  failed  unit  is  being 
repaired.  The  different  availability  states  that  a  parallel 
system  can  be  in  are  shown  Ln  the  following  Markov  graph  in 


?.Z 
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Figure  3.2  (Shccman,  1969:344). 


Figure  3.2  Markov  Graph  for  Two-Unit  Parallel  System 
Availability  with  Repair  (Shooman,  1968:344) 

The  symbols  in  Figure  3.2  are  defined  as  follow*: 

*0  *  both  units  operational 

»  one  unit  failed,  one  unit  operational 

*2  =  both  units  failed 

Xl,x2  ™  unit  designators 

X  =  1/Mean  Time  Between  Failure  cr  1  JO 

//  *  1/Mean  Tim#  To  Repair  or  1/0 

}f  =  2 X  for  an  ordinary  system 
/ 

X  ~  X  for  a  standby  system 

f 

t*  -  /*  for  on*  r«p.iirman 

2  3 


V  *  for  aora  than  one  rtpainas  (k^  >  1) 

•■*  fj  for  ona  repairman 

ft 

U  m  2  fJ  for  two  repairmen 

ft 

(i  *  k2^  for  more  than  two  repairmen  <k2  >  2) 

4  t  *  chang®  in  time 

Mote  that  the  Harkov  graph  accounts  not  only  for  the 
failure  of  zero,  one,  or  both  units,  but  that  it  also 
accounts  for  a  differing  number  of  repairmen  working  on  the 
system.  This  is  an  important  point  to  consider  for  any 
repairable  system.  If  too  many  repairmen  are  working  on  a 
system,  a  point  of  diminishing  returns  is  reached  when  the 
repairmen  begin  interfering  with  each  other's  work.  By 
varying  the  number  of  repairmen  in  the  model,  their  effect  on 
system  availability  can  be  simulated,  and  an  optimum  number 
of  repairmen  can  bo  found. 

Finally,  the  availability  equation  used  in  this  model  is 
derived  from  Figure  3.2  (Shooman,  1968 i 343-346 ) »  The 
equation  is 


K 


t.  • 

I. 


k. 


L 

•  w 

*  » 
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A(t)  «  <1  -  AA/r^)  . 

t  jr  t 

A  A /r^  -  r4(«  3  /r3 

-  'V'V  <i°>  ;■ 
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where 


r3'r4  “  C~<  A  +  X'  + ^  > 

+/-  <  (  A-f-  X  +  /i  +  */  )A 

-  4<  XX  +  X  n'\  /uV)  >**<l/2>  31/2 
and  other  symbol*  arm  as  defined  for  Figure  3.? 


Like  the  availability  of  a 
parallel  system,  the  reliability  is  better  when  the  parallel 
system  has  a  repair  capability.  The  Markov  graph  in  Figure 
3.3  shows  the  different  reliability  states  that  a  repairable 
parallel  system  can  be  in  and  transition  to  (Shooman, 
1968:341). 


Figure  3.3  Markov  Graph  of  Two-Unit  Parallel  System 
Reliability  with  Repair  (Shooman,  1368:341) 
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The  symbols  in  Figure  3c 3  are  defined  as  follows: 
Sq  *  both  units  operational 
s*  **  on®  unit  failed,  one  unit  operational 
&2  m  both  units  failed 

X  =  1/Mean  Time  Between  Failure  or  1/0 

/i  ■  1/Mean  Tim#  To  Repair  or  1/0 
/ 

X  *  2 X  for  an  ordinary  system 


X  “  X  for  a  standby  system 
V  »  /J  for  one  repairman 

iu'  =  k for  more  than  on#  repairman  (k  >  1) 

A  t  *  change  in  time 

Again,  the  Markov  graph  accounts  for  the  failure  of 
zero,  one,  or  both  units  and  the  effects  of  differing  number 
of  repairmen  working  on  the  system.  Note,  however,  that  once 
both  units  are  in  a  failed  state  s  multaneously,  state  s,, 
becomes  an  "absorbing"  state;  the  system  has  reached  the 
limit  of  it's  reliability  for  that  cycle.  Once  the  entire 
system  is  repaired,  the  reliability  cycle  starts  again. 

The  reliability  equation  used  in  this  model  is  also 
derived  from  it's  Markov  graph  (Shooman,  196B : 341-34Z ) . 


R(t)  «  P»Q(t)  +  Ps^t) 


(11 ) 


where 


Ps0<t) 


<  A  t  ft  +  r  /  r  ^ 


r  t 
r  2  1  e  1 


<  A  f  f‘  +  r2f  ~l 


.  r  t 
rv)e  7 


Zb 


P»1(t)  -  {  \'/xl  -  r2)«r1t 

-  <  ^,/r1  ~  rz)9,r ^ 

rlfr2  =  [-<  A  +  A%  u)  +/- 

( (  A+  A'+  U  ')Z  -  4  AA'  >**l/2]l/2 

and  other  symbols  are  as  defined  in  Figure  3.3 

Model .  The  simulation  model  used  the  following  baseline 
parameters : 

Mean  Time  Between  Failure  (1/A  )  =  100  hours 

Mean  Time  To  Repair  (1/ /i  )  *  20  hours 

Simulated  Test  Time  *  175  hours 

Both  the  “true“  sample  MTBFs  and  MTTRs  and  the  simulated 
sample  MTBFs  and  MTTRs  were  drawn  from  the  chi-square 
distribution  in  the  same  manner  described  for  the  single-unit 
system  availability  simulation.  Also,  when  the  same  seeds 
for  the  random  sample  generation  are  used  for  both  the 
availability  and  reliability  simulation,  they  correspond  to 
the  same  sample  systems. 

In  summary,  simulated  availabilities  and  reliabilities 
were  found  for  Type  IT  censored  sample  sizes  of  10,  20,  30, 
and  50.  Point  estimates  were  found  for  the  lower  confidence 
limits  of  0.98,  0.95,  0.90,  0.85,  and  0.B0.  The  model  ran 
for  500  repetitions  with  500  trials  per  repetition. 
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IV .  Monte  Carlo  Simulation  Results 

After  the  Mont®  Carlo  simulation  modal  was  run  as 
described  in  Chapter  3,  the  simulated  point  estimates  were 
compared  against  th*  "true"  lower  confidence  limit  point 
estimates  for  the  single-unit  system  availability  and  the 
two-unit  parallel  system  availability  and  reliability.  The 
success  of  using  this  Monte  Carlo  method  is  based  upon  the 
simulation  coverages  of  the  true  lower  confidence  limit  point 
estimates . 

Single-Unit  System  Availability 

The  Monte  Carlo  simulation  model  was  run  to  find 
estimated  single-unit  system  steady-state  availabilities 
using  the  following  parameters i 

Mean  Time  Between  Failure  =*  100  hours 
Mean  Time  To  Repair  =*  Z0  hours 

Thus,  the  exact  availability  is  as  follows: 

Ass(t)  *  100/(100  +  Z0)  =*  .8333 

The  Monte  Carlo  simulation  coverages  of  this  sing.le-unib 
system  steady-state  availability  are  shown  in  Table  I.  Again, 
all  simulations  were  run  for  500  repetitions  with  500  trials 
par  repetition. 
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Table  I. 

Results  of  the  Monte  Carlo  Method  of  Simulating 
Single-Unit  System  Availabilities 


1  Sample  Size 

T 

Actual 

Coverage  (X) 

1 

i 

.98 

1  .93 

.90 

I  .85 

1 

.80 

1  10 

i 

27.6 

1  93.8 

90.2 

1  86.2 

1 

80. 

1  20 

i 

97.4 

1  95.2 

90.0 

i  86.0 

1 

81.i 

1  30 

i 

96.6 

1  93.4 

88.6 

1  83.4 

1 

7S,i 

1  50 

i 

98.6 

1  95.6 

,  91.2 

I  86.6 

1 

81. 

v« 


Two-Unit  Parallel  System  Availability 

The  Monte  Carlo  simulation  of  two-unit  parallel  system 
availabilities  used  these  parameters: 

Mean  Time  Between  Failure  *  100  hours 
Mean  Time  To  Repair  =  20  hours 

Testing  Time  =  175  hours 

The  exact  availability  (A(t>  where  t  is  175  hours)  for 
the  two-unit  parallel  system  is  0.9460.  (The  exact 
availability  was  found  by  solving  Equation  10  in  Chapter  3 
using  the  above  parameters.)  The  Monte  Carlo  simulation 
coverages  for  this  two-unit  parallel  system  availability  are 
sli  wn  in  Table  II.  All  simulations  were  run  for  500 
repetitions  with  500  trials  per  repetition. 
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Table  II. 


Results  of  the  Monte  Carlo  Method  of  Simulating 
Two-Unit  Parallel  System  Availabilities 


I  Sample  Size  I  Actual  Coverage  <H> 


1 

|  .90 

!  .95 

T 

.90 

T 

.  85 

T 

.80 

j  ib 

i  96.2 

I  93.4 

1 

90.2 

i 

86.0 

i 

81.1 

1  20 

1  98.8 

!  96.4 

I 

92.4 

i 

87.8 

I 

83.  - 

1  30 

I  97.4 

i  94.6 

1 

91 . 4 

i 

86  o  4 

1 

80.1 

1  50 

I  98.4 

1  95.6 

1 

90.4 

i 

85.8 

1 

79.1 

Two-Unit  Parallel  System  Ral ■i.abilitv 

The  Monte  Carlo  simulation  of  two-unit  parallel  system 
reliabilities  used  the  following  parameters! 

Mean  Time  Between  Failure  =  100  hours 
Mean  Time  To  Repair  =  20  hours 

Testing  Time  *  175  hours 

The  exact  reliability  (R(t)  where  t  is  175  hours)  is 
0.6565.  {The  exact  reliability  was  found  by  solving  Equation 
11  in  Chapter  3  using  the  above  parameters.)  Table  111 
presents  the  Monte  Carlo  simulation  coverages  for  the 
two-unit  parallel  system  reliability. 
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Table  III. 


Resu  3  of  t^a  Monte  Carlo  Method  of  Simulating 
Two-Unit  Parallel  System  Reliabilities 


1  Sample  Size  !  Actual  Coverage  ( X ) 


! 

1  .98  1 

.95 

.90 

.85  i 

o 

eo 

1  10 

1  99.6  l 

96.6  1 

92.4 

89.2  1 

1  84. 

1  20 

!  99.4  | 

97.8  | 

94.8 

90.4  | 

I  B7.i 

1  30 

1  98.8  | 

96.2  1 

92.4 

88.6  | 

!  83.1 

1  50 

1  98.8  i 

94.8  | 

89.6 

85.4  | 

1  80.1 

Summary 

In  summary,  Moute  Carlo  coverages  were  obtained  for  the 
lower  confidence  limits  for  single- unit  system  availability, 
the  two-unit  parallel  system  availability,  and  the  two-unit 
parallel  system  reliabi.*  ity . 
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The  primary  abjective  of  this  research  is  to  develc  p  an 
efficient  Monte  Carlo  simulation  method  to  .ind  lower 
confidence  limits  for  the  availability  and  reliability  of 
maintained  systems.  After  analysing  he  results  of  the 
single-unit  system  availability  simulation  and  the  two-unit 
parallel  system  availability  and  reliability  simulation,  the 
conclusion  i*  that  this  Monte  Cax'o  simulation  method  iti  i 
viable  and  efficient  method  of  simulating  lower  confidence 
limit  system  availability  and  readability  point  estimates. 
However,  it  should  be  noted  that  the  coverage  of  system 
reliability  lower  confidence  limits  is  somewhat  high, 
although  still  fairly  good. 

Recommendation*  for  future  research  include  to  further 
study  the  Monte  Carlo  method  of  simulating  ay ft tew  reliability 
lower  confidence  limits,  to  further  study  using  this  Mont* 
Carlo  method  with  varying  parameters  and  different 
situations,  and,  finally,  to  design  and  develop  an 
interactive,  user-friendly  program  which  allows  varying 
input*  for  different  maintained  ay stems. 

In  conclusion,  the  more  rasaarch  performed  on  accurately 
and  efficiently  simulating  lower  confidence  limits  for  system 
availability  and  reliability,  the  more  engineers  can  use  this 
simulation  tool  to  design  and  develop  the  most  effective  and 
reliable  systems  and  equipment  possiL  & .  Uit  imately ,  the 
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batbar  aquipmant  and  aystam*  fialdad  to  U,  S.  Armed  Force® 
ptnoantl  F  th«  hatter  bha  country '  ss  defensive  and 
war-fighting  capabilities  will  b«  now  and  in  the  future. 


Appendix  A 

gjiBB^laa  MiAhfld 

As  mentioned  in  Chapter  3,  "The  Montu  Carlo  Simulation 
Modal , *  sample  0s,  or  Mean  Tims  Between  Failures  (MTBFs) ,  and 
sample  <£»,  or  Mean  Time  To  Repairs  (MTTRs) ,  were  drawn  using 
the  chi-square  distribution  as  follows: 

©  -  (x22r0)/2r  <A1' 

where 

2 

X  *  random  number  drawn  from  the  chi-square 
distribution  with  2r  degrees  of  freedom 
r  «  sample  size 

MTTRs  were  generated  using  the  same  equation  with  sub¬ 
stituted  for  0. 

Ma».y  methods  of  generating  random  numbers  from  the 
chi-square  distribution  are  available,  but  only  two  methods 
were  used  in  this  research.  The  first  method  was  to  directly 
generate  chi-square  distributed  random  numbers  using  the 
International  Mathematical  Statistical  Library  (IMSL)  GGCHS 
(chi-square  random  deviate  generator)  subroutine.  The  second 
method  vised  was  to  generate  a  random  number  using  the 
chi-square  distribution's  r«productiv«  property  (Miyamura, 
1902:315).  This  method  involved  summing  tha  absolute  value 
cf  the  natural  logarithm  of  a  uniform  (0,1)  random  number  r 
(r  is  the  sample  size)  times.  The  uniform  random  number  was 


generated  by  the  IMSL  function  GGUBFS  which  is  a  basic 
uniform  (0,1)  random  number  generator  function. 

There  is  virtually  no  difference  between  the  two  methods 
of  generating  chi-square  distributed  random  numbers.  Two 
simulations  of  sample  size  10  were  run,  each  using  a 
different  generation  method  and  the  same  seed.  The  results 
are  shown  in  Table  1A  and  Table  IIA. 

Table  IA. 

Comparison  of  Coverages  of  the  Two  Random  Number 
Generator  Methods  for  the  Single-Unit  System  Availability 


1 

Method 

!  Actual 

Numbers  and 

Coverage 

“> 

1 

1  .98  | 

.95  i 

.90  1 

.85  | 

.80 

1 

GGCHS 

I  97.6  | 

93.8  i 

90.2  | 

86.2  1 

80.0 

I 

I  488  1 

469  | 

451  1 

431  | 

400 

1 

GGUBFS 

I  97.6  | 

93.8  S 

90.2  | 

86.2  | 

GO.  0 

1 

1  488  | 

469  | 

451  | 

431  | 

400 

Table  I I A 
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Comparison  of  Coverage*  of  the  Two  Random  Number  Generator 

i 

Method*  for  the  Two-Unit  Parallel  System  Availability 


! 

Method 

1  Actual 

Numbers  and 

Coverage 

”"<*>  i 

p  "T* 

r 

1 

1  .98  i 

.95 

.90 

.85  j 

.80  1 

! 

GGCHS 

1  96.2  t 

93.4 

90.2 

86.0  1 

81.6  | 

1 

1  481  1 

467 

451 

431  I 

408  | 

1 

GGUBFS 

1  96.2  { 

93.4 

90.2 

86.0  1 

81.8  1 

s* 

1 

1  481  1 

467 

451 

431  1 

409  1 
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PROGRAM  AAAAA 


INTEGER  A,B,C,D,N,E,R 
REAL  TTIME , MTBF , MTTR  ,  CHI  2 

REAL  LAMBDA, MU 

REAL  AROTS(3QO,3,3) 

REAL  ANTER { 300 , 3  > 

REAL  LAMBPR 

DOUBLE  PRECISION  DSSID 

DATA  AR0TS/750Q*0.Q/ 

DATA  AMTSR/1300*0 . 0/ 

OPEN (UNIT*  , FILE*'  ') 

REWIND (UNIT-  ) 

DSEED  *  XXXXXXX 
TTIME  *173 
N  *  XX 
R  -  XX 

DO  1G0  A  -  1,300 
CHI 2  -  0.0 
DO  300  E  -  1,R 

CHI 2  -  - <LOG(GGUBFS( DSEED) > >  +  CHI 2 
300  CONTINUE 


CALL  GGCHS ( DSEED , N , R , CHI 2 ) 

MTBF  *  <CHI2*100)/N 
DO  600  E  -  1 ,  R 

CHI 2  *  -( LOG (GGUBFS( DSEED) ) )  +  CHI 2 
600  CONTINUE 

MTTR  -  (CHI 2* 20  )  /N 
DO  200  B  -  1,300 

CALL  GNEXP ( DSEED , LAMBDA , MU , MTBF , MTTR ) 

LAMBPR  *  LAMBDA  A  2 

B2  *  -(LAMBDA  +  LAMBPR  ♦  MU  +  MU ) 

AC1  »  ( ( - ( B2 ) )**2> 

AC2  *  <1*  <  LAM H DA* LAMBPR  +  LAMPRAMU  +  MU*MU  ) 
AC  -  SQRT ( AC1  -  AC2) 

R3  *  ( B2  +  AC ) / 2 


PARI  »  <l-< <LAMBDAALAMBPR)/(R3AR4) ) ) 

PAR 2  *  ( LAMBDAALAMBPR } / ( R3-34 ) 

PAK3  -  <  (EXP (R3ATTIMS)  )/R3)~<  <H3<P<R4ATTIME )  )/R4) 
AVAL  *  PARI  -  <  PAR2*PAR3 ) 


LAMBFR  *  LAMBDA  *  2 

B2  -  -(LAMBDA  +  LAMBFR  +  MU) 

AC  -  SQRT< ( <~(B2> > **2 > - < 4*LAMBDA*LAMBPR ) ) 

R1  »  ( 32  +  AC)/2 
R2  -  <B2  -  AO/2 

NUMi  »  ( LAMBDA+MU+R1 ) * ( EXP ( R1 ATTIME ) ) 

NUM2  -  ( LAMBDA+MU+R2 ) *  <  EXP ( R2ATTIME ) ) 

DEN  *  R1  -  R2 

PSO  -  (NUMI /DEN )  -  (NUM2/DEN ) 

NUM3  *  ( (LAMBPR/DEN)A(EXP(R1ATTIME) ) ) 

NUM4  -  ( (LAMBPR/DEN)A(EXP(R2ATTIMK) ) ) 

PS1  *  NUM3  -  NUM4 
REL  »  PSO  +  PS1 

OR 

R5  »  MU/ (LAMBDA  4-  MU) 

Rb  »  LAMBDA/ ( LAMBDA  +  MU) 

AVALS  «  RS 

ANTER(B,1)  -  R1  OR  R3  OR  RS 
ANTER(B,2)  «  R2  OR  R4  OR  R6 
ANTER (  B ,  3  )  »  AVAL  OR  REL  OR  AVALS 

IF  (B  .EQ.  500)  THEN 
CALL  ORDER (ANTER) 

CALL  RROOTS ( ANTER , A , AR0T3 ) 

ENDIF 

200  CONTINUE 

100  CONTINUE 

FORMAT (T2 , F20 . 15 , 2X , F20 . 15 , 2X , F20 „ 15 ) 

DO  300  C  ■  1,500 
DO  400  D  =  1,5 

WRITE (  UNIT*  ,FM'T»  )  AROTS  ( C  ,  D  ,  1  )  ,  AROTS  ( C  ,  D ,  2  )  , 

AROTS ( C , D , 3 ) 

400  CONTINUE 

300  CC..a'INUE 

REWIND  ( UNIT3*  ) 

CLOSE (UNIT-  ) 

STOP 

END 
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SUBROUTINE  GN EXP  <  DSEED , LAMBDA , MU , MTBF , MTTR ) 

INTEGER  N,E,R 
REAL  LAMBDA, MU 
REAL  MTBF, MTTR 
REAL  CHI 2 
REAL  THETA, PH I 
DOUBLE  PRECISION  DSEED 

N  -  XX 
R  »  XX 

CHI 2  »0,0 
DO  700  S  -  1,R 

CHI 2  »  - <LOG<GGUBFS< DSEED) ) )  +  CHI 2 
700  CONTINUE 

THETA  •  <CHI2*MTBF)/N 
LAMBDA  «  1/ THETA 

CHI 2  «  0.0 
DO  BOO  R  «  1,R 

CHI 2  «*  - < LOG ( GGUBFS < DSEED ) ) >  +  CHI2 
BOO  CONTINUE 

PHI  -  (CHI2*MTTR)/N 
MU  -  1/PHI 

•  RETURN 

END 
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SUBROUTINE  ORDER < INTERM) 


INTEGER  C, SWITCH 
REAL  INTERM (500, 3) 

REAL  TEMPI, TEMP 2, TEMP 3 

SWITCH  -  1 

10  IF  (SWITCH  .EQ.  1)  THIN 
SWITCH  »  0 
DO  100  C  *  2 , uOQ 

IF  ( INTERM ( C , 3 )  .GT.  INTERM (C-l , 3 ) )  THEN 
TEMPI  *  INTERM (C,l) 

TEMP2  «  INTERM (C, 2) 

TEMP3  »  IN TERM ( C , 3 ) 

INTERM ( C , 1 )  -  INTERM ( C-l, 1} 

INTERM (C, 2)  *  INTERM ( C-l , 2 ) 

INTERM (C, 3)  »  INTERM (C-l, 3) 

INTERM ( C-l, 1)  «  TEMPI 
INTERM (C-l ,2)  -  TEMP2 
INTERM ( C-l ,2>  «  TEMP3 
SWITCH  -  1 
ENDIF 

100  CONTINUE 

GOTO  10 
ENDIF 

•  • 

RETURN 

END 


4  0 


SUBROUTINE  RROOTS (RINTER , A f ROTS ) 

INTEGER  A,C,D 

REAL  RINTER <300, 3) , ROTS( 500 , 5 , 3 ) 

REAL  LCB , PRQPQR , ROOT1 , ROOT2 , I 

REAL  POLATE , DIFF 

DO  100  D  *  1,3 

IF  <  D  .EQ.  1)  THEN 
I  »  490.632 
ENDIF 

IF  ( D  .EQ.  2)  THEN 
-  473.68 

ENLaF 

IF  (D  .EQ.  3)  THEN 
I  -  430.66 
ENDIF 

IF  (D  .EQ.  4)  THEN 
I  *  425.64 
ENDIF 

IF  (D  .EQ.  3)  THEN 
I  *  400.62 
ENDIF 

DO  200  C  »  1,300 

LCB  -  RINTER <500, 3) 

ROOT1  »  RINTER <300,1) 

ROOT2  «  RINTER <300, 2) 

IF  (I  , LT .  C)  THEN 
PROPOR  -  I  -  <C-1) 

POLATE  *  RINTER (C, 3)  -  RINTER (C-l, 3 
DIFF  -  PROPOR  *  POLATE 
LCB  -  DIFF  +  RINTER (C-l, 3) 

ROOT1  *  RINTER < C-l, 1) 

ROOT2  -  RINTER < C-l ,?,) 

GOTO  300 
ENDIF 

200  CONTINUE 

300  ROTS  ( A  ,  D  ,  1  )  =»  KQOT1 

ROTS ( A , D , Z )  -  ROOT2 
ROTS ( A , D , 3 )  -  LCB 

100  CONTINUE 

RETURN 

END 
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